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SUMMARY 


A new approach to the calculation of the dynamic characac- 
teristics of the combined man-manii ulator-system in manual air- 
craft control has been derived from a model of the neuromus- 
cular system similar to that described by McRuer and Magdaleno. 
(Ref. 1) This model combines the neuromuscular properties of 
man with the physical properties of the manipulator system 
which is introduced as pilot-manipulator model into the manual 
aircraft control. The assumption of man as a qiaasilinear and 
time- invariant control operator adapted to operating states 
- depending on the flight phases - of the control system gives 
rise to interesting solutions of the frequency domain trans- 
fer functions of both the man-manipulator system and the 
closed loop pilot-aircraft control system. It can be shown 
that it is necessary to introduce the complete precision pilot- 
manipulator model into the closed loop pilot-aircraft trans- 
fer function in order to understand the well known handling 
quality criteria of MIL-F-8785B/C , and to derive these criteria 
directly from human operator properties. 


INTRODUCTION 


The pioneer work on the precision pilot model presented 
by D.T. McRuer and h's co-workers (Refs. 1...3) has become the 
most important step towards our understanding of the role of 
man in manual control. It is the combination of neurophysio- 
logy, physical dynamics, and control theory which gives the 
fascinating aspects of how the several complicated problems 
of manual control should be solved. But, there is a little gap 
between the theory and practical solutions, because the pre- 
cision theory always turns out to be too complicated if en- 
gaged to solve such problems as the question for the best mani- 
pulator characteristics in a high performance aircraft. 


In fig. 1 the precision pilot model is shown as it was pre- 
sented by McRuer and Magdaleno (Ref. 2), which enlightens sche- 
matically, what is running in the neuromotor system when the 
pilot puts his (land (or legs) on the manipulator. For practi- 
cal use however, the human engineer wants to revise this model 
to also parametrically based but simpler facts. 
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Fig. 1: Precisiion pilot model, after McRuer and Magdaleno 
(ref. 1) 

If you ask the pilot about manual control qualities, he ex- 
plains his wishes on 

stick forces 

aircraft response (as it is) 

predictability of response (as he expects it is) 
lead elements (after inquiry) , 

but he is not accustomed to consider his stick force and dis- 
placement feedback. If the pilot must consider these, the 
handling qualities of the system may be bad. 

So for practical use we tried to reorganize the pilot model. 
Fig. 2 shows another pilot model less sophisticated as that 
of fig. 1 but more practical in use. Force and displacement 
feedback now feeds into the spinal chord without becoming con- 
scious to the operator who is engaged in the compensation of 
system lags by lead. This model too can be fully identified 
by physical parameters. (Fig . 3) We now can divide the pilot 
model again into two parts; part A identifies the mental para- 
meters Uead) while part B represents the neuromuscular-physi- 
cal parameters of the combined man-manipulator system (lag) . 
Now the intention of the following analysis is the formulation 
of the frequency domain transfer function of this quasilinear 
pilot model and the calculation of amplitude and phase shift 
of the nan-manipulator combination? while further important 
investigations include the pilot-aircraft closed loop charac- 
teristics and their influence on handling qualities. (Ref. 5) 
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ANALYSIS 


Formal arrangement of the model. 


The model (fig. 2) can be arranged formally alike 


( 1 ) 

F ( S ) 

F^(s) .F^(s) . Fj_, (s) .F 



P 

1-tF (s) ,F, (3)+F„(s) 
lx x> w 

, where 

(1.1) 

F ( S) = 

A 

- c s 

K . 1 +T. s ) e A 

A A 

(information input term) 

' 1 . 2 ) 

F^(S) - 

_ r yS 

K . ( 1 +T f ) e 
V V 

(information processing term) 

(1.3) 

F^(s) . 

•^E 

UR^S.H^S' 

(second order lag of man- 
manipulator) 

(1.4) 

F (b) = 

Li 

e 

(neuronal impulse delay) 

( 1 - S ) 

Fj^(s) = 

( 1+T|^s)e 

(delayer lorce-feedback ) 

(1.6) 

= 

_Cw^ 

(delayed displacement feedback) 


Rearranging the equations (1 . 1 ) . . . (1 .6) into eq. (1) results 
in 
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and rewritten by the fol levying introductions 
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To achieve the roots of this system, we first have to ell- 
minate the delay terms as much as possible > which is done by^ 
multiplication of both the numerator and denominator of eq. (^) 

by This results in the replacement of e ^ by e 

e ^ ® ® in the numerator/ and another form of the denominator 

D <s) of eq. (3) like 
P 


(4) D (s) = ( 

P 


1+RgS + MpS’ )'[ e''^®+Kj^(l+Tj^s)]+KgK^(l+T^s) 


While the numerator has already its final root arrangement, 
the denominator has not. For better Handling/ the term 
s 

e should be approximated by a Taylor series 

(5) = 1 S+C' +R ( s ) 


2 ! 


(5.1) R (s) = 


n = 2 , 3 , 6 


the frequency range of which is 0< (s = j )< 70 rad/sec. with 
an amplitude and phase error below 20 %/ for n = 2: 


(6) F (.s) = A(s)B(s) = 
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K ( 1 t-T s ) ( 1 +T, .s ) e 
p A V 


TeS 


( 1 t-R -S + M sM fl + f ( 1 +Tj^r, )J (1+T^s) 
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The roots of the denominator polynom are solved by the calcu- 
lation of the denominator polynominal form 

(7) D = a„+fi.s + .j,,s''+.is’+d.s'l+ara^ 
p i 3 4 J 

and by comparison of the coefficients a. and the coefficients 
of another/ well-known polynom the roots of which meet the 
5th order arrangement 

(l+b,s) (l4-b,s+b.s'‘ ) (1+b s+b s' ) 

1 2 i 4 5 

[h" necessary arithmetics are described and discussed in ref. 5 
It has to be mentioned/ that an exact solution of a 5th order 
linear arithmetic equation does not exist. 


SOLUTION 


The best solution of eq. (6) is the following 
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which describes the required roots of a 5th order lag system 
B(s), if 

-r G 

( 9 . 1 ) A(s) = K ( 1 +T s) ( 1 +T. s)e ® 
p A V 

The conditions for the solution of B(s) within eq. (9) are; 
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is a "high frequency weighting factor", in which , k, and n 
are defined as follows 
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n = 2 (see eq, 5.1.) 
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The 5th order .solution (9) degrades to a 3rd order solution, 
if 6'>0. A limit of S is given by&<1: If 5 <£. set 0 , whiles- 
is defined by 


(10.9) 


i . e . if 




q;’ = 1 > 70 rad/sec which is the limit in u» de- 


/kii-M 


fined by the approximation in eq. (5.1) 



Numerical evaluation 


While the numerical values of Mjj, (limb inertia) , Mg (mani- 
pulator inertia) , Cg (feel spring constant) , and Rg (manipula- 
tor viscous resistance) can be achieved by measurement - or 
effectively are Icnown - the values of Cm (limb muscular tension 
at operating point) or Rm (limb viscous resistance at operating 
point) have to be assumed in order to achieve reasonable results. 


The same is done for the values of Kj^ and Ky, which only 
can be defined statically as in the following way 

K, = (P = control sticlc force) 

^ F max ® 

s 

K = 6s point) _ control sticlc displacement) 

6-g max ® 


while Kp might be defined as 


Kg = max 


(F /n = stick force gradient at 
operating point, eg. n = n for 
pitch axis) ^ 


For force-displacement-manipulators, the product K K is de- 
fined to be unity at every operating point. 


RESULTS 


Comparison with experimental data 


Some calculations of T^^, and are made based on 

experimental data publishea in ref. 2. Most interesting data 
were those of ref. 2 which have been evaluated from tracking 
experiments with all manipulator characteristics but only the 
simplest controlled element (pitch axis) 

Fc(s) = K^, 

because these data are most characteristic for the man- manipu- 
lator system alone. If one assumes K^=Ky=1 and T^ -Ty-#' 0 for a 
controlled element without the necessity of lead compensation, 
which will be almost true for open loop control, the resulting 
transfer function is 

F(s) = B(s) .Kc, Kc = 1 . 

Using the data of experimental runs no. 3, 8, and 23 of 

ref. 2 for M , R , C and assuming: 
s s s ^ 
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run no. J 



Fig. 4: Comparison between experimental data (®) measured by 
Magdaleno and McRuer (ref. 2) and theoretical data (-) 
- pitch axis, open loop - 





! 


M =0.0.0199 Tkam^J = hand-to elbow inertia for side stick 
m 

C =1.5 rNmVsec*J 

m L / j 

(Mm+Ms) (Cm-*-Cs) (see ref. 5) for X =0.7 

’"'2Ssf^:nr 

the calculated open loop Bode characteristics are plotted 
against the experimental data from ref. 2 in fig. 4. As one 
can recognize, only the amplitude differences are quite 
noticeable, because the calculation was based upon ^nd 

Kp = 0,666 for run 3 (force-displacement stick) 

0,666 •' " 8 (displacement stick) 

0,00913 " " 23 (force stick) 

Nevertheless, the calculated plots in fig. 4 clearly show 
the same characteristics as the experimentally evaluated ones. 


Handling quality requirements 


According to the Handling quality specification MIL-F-8785-C 
(ref. 6) the phase angle between Fg (stick force) and <5“^ (con- 
trol area displacement) should not exceed-35 degrees, for 
level 1, related to the axis resonant frequency , £J n* The mini- 
mum requirements for 1 /Tw and 5 e^E ^re, while other phase 
shift between Fs and 5 be excluded. 


= (3, 5. ..7) for ? 1.0 

Assume, the pitch axis resonant frequency 6/ = 3... 5 

rad/sec, these minimum requirements mean: 

^ = (3. ..5). 3, 5 = 10. ..17,5 rad/sec (T = 0,057. ..0,1) 

U; = ( 3 . . . 5) . ( 3 , 5 . . . 7 ) = 10. ..35 rad/sec 

As an example, we should remember the high resonant fre- 
quency of the T33 manipulator system, used by CAL during the 
Neal/ Smith experiments, which was it/ = 31 rad/sec (ref. 7). 
The low values of as well as the^high values of u» g from 
the above assessment are realized only by light, stiff mani- 
pulators (hand controls) used aboard fighter aircraft with 
high resonant frequencies ((t^nsp ~ 3... 10 rad/sec). 







-A15- 



Handling quality criteria derivations 


The derivation of well known handling quality criteria is 
possible by the introduction of the solution of eq. (9) into 
the closed-loop transfer function of the system pilot-air- 
craft. With the assumption thatii^E is high enough to satisfy 
the requirement 


V E - 3 , 5 . . . 7 ) .u; 


nsp 


- pitch axis - 


the solution of eq. (9) is simplified to 
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If this pilot model is used in a closed loop together with the 
well known pitch axis aircraft transfer function 
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the closed loop response is 
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is the best substitution. The solution of the closed loop 
problem (13) can be achieved by using a method similar to that 
applied to eq. (6) . We may suggest 
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These conditions show a very strong dependence between the 
aircraft parameters T » I) the most critical para- 

meter of the man-manipulator T^r and the lead time constants 
generated by the pilot himself in order to stabilize the closed 
loop. If level 1 proven aircraft parameters are used, the 
closed loop response satisfies the Neal and Smith tracking 
criterion (ref. 7) for good handling, if the Bandwidth of 
the "band-pass filter" described by eq. (15) is restricted to 

BW = 3,5 rad/sec 
(BW) = - 3dB 

The pilot conditions are then: 

T <0,4 sec (PlO-criterion) , k v. Fs 
® ' ni 

T +T^ < 2,0 sec (ref. 8) 

V A 

This "band-pass filter criterion" is shown by fig. 5. While 
Tv is the indicator for critical a/c parameters, Ta is that for 
critical man-manipulator parameters. Both should be as small 
as possible for good handling qualities. 

The addition of the 2nd order term ( 1 +kRES+kMj;S* ) disturbes 
this Bandwidth conditions only, if VvHcMe « < 4 »nsp* In this 
case, another lead time constant TaI has to be generated by 
the pilot in order to compensate the lag time constant kRg, 
resulting in a higher order pilot lead element 


(UT^s) (l+T^^s) (1+T^s) , 

the sum of which again must satisfy 


T^+T^^+T^< 2.0 (Arnold, ref. 8). 

The higher order of the lead element surely will also degrade 
the pilot rating further. 


Fig. 5: Handling 
quality criterion pro- 
posed by Neal and 
Smitli (ref. 7). Bode 
diagram of closed 
loop response of 
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CONCLUSIONS 


The formal arrangement of the full quasilinear pilot model 
first proposed by McRuer et. al . into both open loop and 
closed loop transfer functions of the manual aircraft control 
is a powerful means to assess aircraft handling quality cri- 
teria. The derivation of known criteria from lead term and de- 
lay term limits of the pilot has been shown clearly. 
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